Single photons from relativistic collision of lead nuclei at CERN SPS: A reanalysis 



o 
o 

(N 



> 

cn 
o 

O 

o\ 
o 



X 



Rupa Chatterjee,^ Dinesh K. Srivastava,^ and Sangyong Jeon^ 

' Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700 064, India 
''Department of Physics, McGill University, Montreal, Canada H3A 2T8 

We present a reanalysis of single photon production from relativistic collision of lead nuclei at 
CERN SPS measured by the WA98 experiment. The refinements include use of iso-spin and shad- 
owing corrected NLO pQCD treatment for prompt photon production using an optimized scale for 
factorization, renormalization, and fragmentation and use of hydrodynamics suited for non-central 
collisions along with a well tested equation of state admitting a quark-hadron phase transition. A 
quantitative explanation of the data requires a large initial temperature (at a small formation time 
of about 0.2 fm/c) and a moderate increase in the prompt yield which could perhaps be attributed 
to the Cronin effect in nuclei. The data can also be explained using a moderate initial temperature 
(at a formation time of about 1 fm/c) with a very large ii'-factor multiplying the prompt yield. We 
show that different initial times give rise to different values for the elliptic flow parameter 172 for 
thermal photons. We also show that a measurement of V2 for thermal photons could also distinguish 
between the scenarios with or without a phase transition. 



I. INTRODUCTION 

The first observation of single photons in relativistic 
collision of lead nuclei in the WA98 experiment at CERN 
SPS [Ij remains an important mile-stone in our search 
for the quark-hadron phase transition. The earlier ex- 
periment studying the S -\- Au collisions had provided 
only (though quite useful) upper limit on the single pho- 
ton production ^ . The importance of the single photons 
stems from the expectation that once produced they leave 
the system without any further interaction (see Ref 3, 
for a recent account of nodal developments in this field). 
It is thus, expected that if a thermalized system of quarks 
and giuons or hot hadrons is produced in such collisions, 
its temperature could be related to the spectrum of the 
single photons. On the experimental front, the success of 
this endeavour hinges on our ability to subtract out the 
decay photons from the inclusive spectrum of photons, 
while on the theoretical front it depends on our ability to 
evaluate non-thermal photons in a quantitative manner. 

In the present work, we re-analyze the single photon 
measurements reported by the WA98 experiment after 
incorporating several recent improvements in our under- 
standing of the physics of relativistic heavy ion collisions. 
Firstly, we perform the NLO pQCD evaluation of prompt 
photon production using the optimized scale for factor- 
ization, renormalization, and fragmentation, Q — pt/2 
which has been found to describe a vast body of sin- 
gle photon production in pp collisions without introduc- 
tion of any intrinsic fcy for protons [1, [H . We explicitly 
account for the iso-spin of the projectile target system, 
which affects the results at large xt = '2.pTl\fs and in- 
clude the effects of parton shadowing @. Next we ac- 
count for the azimuthal anisotropy of the system for non- 
central collisions, while performing the hydrodynamics 
calculations. Finally we explore the effects of varying 
the initial conditions to set limits on the likely initial 
temperature. 

In the next section we briefly discuss our estimates for 
the prompt photon production. In Sect. Ill we describe 
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FIG. I: Prompt photons from pp collisions at ^/s = 19.4 GeV. 
Experimental results for E704 7] for pp collisions and those 
estimated from p-|-^^ C collisions by the E629 j| and NA3 Q 
experiments are also given for a comparison. The results for 
the NA3 experiment use two different triggers; conversion and 
calorimeter. 



the setting up of the initial conditions, and in Sect. IV 
we discuss the results. Finally we give our conclusions. 



II. PROMPT PHOTONS 

As mentioned earlier, it is quite crucial to get an ac- 
curate estimate of prompt photon production in nucleus- 
nucleus collisions, before we can begin to explore the ini- 
tial conditions of the thermalized system. The success 
of PHENIX experiment in measuring the single photon 
production at large pt in Au -I- Au collisions at RHIC 
energies has brought this consideration into a sharp fo- 
cus. Thus, for example, it is now realized [lo| that the 
'suppression' of single photons at large in Au-V Au col- 
lisions compared to those from pp collisions at the same 
nucleon-nucleon centre of mass energy has its origin pre- 
dominantly in the difference of isospin for protons and 



10 fTT 



A+A, 17.3 AGeV 




1.00 F 



0.95- 
0.90 - 
J 0.85 - 
O80 - 



10 



1.2- 



sa 0.8-^ 



0.6- 



10" 



3.5 4.0 4.5 
Pj (GeV/c) 



— Pb+Pb, 5500 AGeV 

— Au+Au, 200 AGeV 

— Au+Au, 62.4 AGeV 

— Pb+Pb. 17.4 AGeV 




Prompt y, NLO pQCD 
Lso-spin 



10" 



10 



10 



x^=2p^s 



— Pb+Pb, 5500 AGeV 

— Au+Au, 200 AGeV 

— Au+Au, 62.4 AGeV 
Pb+Pb, 17.4 AGeV 



Prompt y, NLO pQCD 
Iso-spin & Shadowing 




10' 



10 



10 




Pb+Pb @SPS 
b = Ofm 



10 



0.4 



8 10 
T (fm/c) 



12 



14 



16 



0.3 



o 

A 
H 
V 



0.2 

0.1 



0.7 



0.5- 

A 

^ 0.4 

H 
> 

V 0.3 



0.2 
0.1 
0.0 



8 10 
T (fm/c) 



12 



14 16 



6 8 10 
T (fm/c) 



12 14 



16 



J 1 1 M M 1 M M 1 M 1 1 M M 1 1 M 1 1 M 1 M 

I Pb+Pb @SPS 


1 1 1 M 1 1 M 1 M M 1 1 M 1 M M 1 M M 1 1 [ 


1 b=0 fm 


<«T„= 0.2 fm/c \ 




Tg= 0.4 fm/c : 




^ T„= 0.4 fm/c -1 




^T„= 0.8 fm/c - 




„ Tg= 1.0 fm/c : 


T, 1 1 1 In 







"1 1 




; Pb+Pb @SPS 






b=0 fm 








^^^^ "V 


0.2 fm/c : 






0.4 fm/c : 






0.6 fm/c = 






0.8 fm/c : 






1.0 fm/c J 




M.I 1 





x.j,=2p^s 



FIG. 3: Evolution of average energy density (upper panel), 
temperature (middle panel), and radial flow velocity (lower 
panel) with time for different initial times to but identical 
rapidity density for total entropy and net baxyons for a central 
collision of two lead nuclei at SPS energies. 



FIG. 2: Upper panel: Effect of iso-spin and parton shad- 
owing on production of prompt photons, calculated using 
NLO pQCD, at ^/sNN ~ 17.3 GeV, which corresponds to 
the nucleon-nucleon centre of mass energy for the WA98 0] 
experiment. Results are given in terms of the nuclear modi- 
fication factor Raa for pn, nn and PbPh collisions. Middle 
panel: Effect of iso-spin at SPS, RHIC, and LHC energies as 
a function of xt ~ 2pT / y/s. Lower panel: Same as before 
with shadowing, (see Ref.[ll|) 



neutrons (or their valence quark structure). 

This is often overlooked when the pp data are scaled 
by the nuclear thickness TaaI^), for the above com- 
parison. Of course, one additionally needs to account 
for effect of jet energy loss if a quark-gluon plasma is 
formed. The study of prompt photon production in pp 
collisions has reached a high degree of sophistication. All 
the available data have now been analyzed using NLO 
pQCD H i, [H, m and it is generally believed that 
choosing the factorization, renormalization, and frag- 
mentation scales as equal to pt/2 provides an excellent 
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FIG. 4: Evolution of average effective temperature with time 
for different initial times to but identical rapidity density for 
total entropy and net baryons for a central collision of two 
lead nuclei at SPS energies. 



description to all the single photon data except for those 
from the E704 and the E706 experiments, without 
the requirement of any intrinsic kx ■ Inclusion of intrinsic 
kx improves the description of the these data but simul- 
taneously destroys the good agreement with all the other 
data. The E704 data is at ^ = 19.4 GeV, which is close 
to the y^SNN ~ 17.3 GeV relevant for the WA98 experi- 
ment. Two other experiments NA3 9] and E629 8] have 
measured single photon production from p C colli- 
sions at the same energy and these data are often used 
with a normalization by the mass-number of the target to 
estimate the pp data, even though half of the nucleons in 
the target are neutrons. We have verified that accounting 
for this reduces the theory values by about 2% at pT ~ 
2 GeV/c and by about 15% for p-r ~ 6 GeV/c, which 
is well below the other experimental uncertainties. We 
shall ignore this for the moment. We show our calcula- 
tion for the single photon production at •\/s= 19.4 GeV in 
Fig. [2 along with the data reported by the NA3 [1], 0], 
and E629 [8| experiments. We note that the fragmenta- 
tion contribution at this energy is of the order of 30% 
of the Compton -I- annihilation term. We also see that 
the NLO pQCD provides a good description of the NA3 
data, while it underestimates the E704 and E629 data by 
a factor of 2-6. This has been noted by several studies as 
mentioned earlier [l^ and it is known that these data 
deviate also from the xt scaling which all the other data 
follow |15|. We have verified that this scaling is in good 



agreement with the NLO pQCD results for values of px 
up to about 4.5 GeV/c, but over-predicts the results con- 
siderablyat higher pT- We are discussing this point again 
as several studies have tried to accommodate these data 
by incorporating intrinsic kx for the partons, which is not 
favoured by the rest of the data. We must add though 
these results are among the earliest measurements of sin- 
gle photons, which may account for the inconsistency of 
data between different experiments and even within the 
same experiment ^9,]. 
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FIG. 5: Upper Panel: Fit to single photon spectra from 
Pb(158 AGeV)+Pb collisions measured by the WA98 1] ex- 
periment for To = 0.2 fm/c The prompt photon contribution, 
is scaled by a factor of 2.7 to normalize the theoretical re- 
sults to the experimental data at pt = 2.55 GeV/c. Middle 
panel: Details. Lower panel: Elliptical flow coefficients for 
the thermal photons. QM and HM stand for photons from 
quark matter and hadronic matter. 



Recently, prompt photon production in p C and 
p_l^208p^ collisions at CERN SPS energy appropriate for 
the WA98 experiment [T^ has been measured. Only the 
upper limit of the single photon production could be de- 
duced. We have verified that the upper limits are about 
a factor of 5-10 larger than the NLO pQCD calculations, 
though the slope of the data is described well by the cal- 
culations. 

We also note that the inclusion of intrinsic kx of par- 
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FIG. 6: Same as Fig.[5]for tq = 0.4 fm/c. Note that the NLO 
pQCD results have to be scaled up by 4.8 for describing the 
data 



FIG. 7: Same as Fig.[5]for to = 0.6 fm/c. Note that the NLO 
pQCD results have to be scaled up by 5.4 for describing the 
data 



tons is not easy even at lowest order of pQCD [T3, 
and the results for NLO pQCD are often inferred by us- 
ing a K-factor which describes the difference of results 
of lowest order pQCD with and without the intrinsic fc^ 
(see also Ref. ^]). In a nuclear medium, Cronin effect 
also contributes to the broadening of the transverse mo- 
mentum spectrum. In Ref. [l^l , it was reported that this 
broadening by the Cronin effect can lead to an enhance- 
ment of photon production by a factor of about 2.5 for 
Pt of 2-4 GcV/c for the case of the WA98 experiment. 

Our task of obtaining an accurate estimate of prompt 
photon production in Pb + Pb collisions at y^ijvTV — 17.3 
GeV is further complicated. Firstly, there is no single 



photon production data from pp collisions at this energy. 
In any case there is no data for pn and nn collisions, 
which will also contribute to the production of prompt 
photons from lead nuclei. The importance of these can 
be seen from Fig. [5] where we show our results for the 
effect of isospin and shadowing for prompt photon pro- 
duction at this energy. We see that effects of isospin and 
shadowing reduce the single photon production at pr = 4 
GeV/c by about 30% for Pb + Pb collisions compared to 
naive scaling of pp data by the nuclear thickness used in 
most of the early studies, including those involving one of 
the present authors [llj . We also note that the inclusion 
of shadowing leads to a significant variation in the end 
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FIG. 8: Same as Fig.[5]for tq = 0.8 fm/c. Note that the NLO 
pQCD results have to be scaled up by 5.7 for describing the 
data 
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result so that at lower transverse momenta, the single 
photon production goes up. The middle and the lower 
panels of the figure show the results for xt = 2pt/\/s 
scaling with the inclusion of iso-spin and parton shadow- 
ing at SPS, RHIC, and LHC energies[ll|. The deviations 
from a scaling behaviour are due to the scale dependence 
of the structure functions and the QCD interactions [2^ . 

In the light of the discussions above we take the fol- 
lowing view for getting the yield of prompt photons for 
Pb + Pb collisions at an energy corresponding to the 
WA98 experiment. For a given impact parameter 5, we 
first estimate the effective number of protons and neu- 



trons from the number of participants: 



where 



Npi.-rtib) ^ / dxdyv{x,y,b) 



'^{x,y,b) = 
[TAix+b/2,y) 

+TB{x~b/2,y) 



(1) 



l-{l~aTBix-b/2,y)/By 



{l-aTA{x+b/2,y)/Af]}{2) 



is the surface density. In the above Ta is the nuclear 
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FIG. 10: Upper Panel: V2 for thermal photons for different 
To. The resuhs for photons Irom hadronic matter alone are 
also given. Lower Panel: v-^ for single photons. Results below 
Pt = 2.55 GeV/c (shown as dotted curve) are obtained by 
arbitrarily using the normalizing lactor at pr = 2.55 GeV/c. 
It is assumed that the difference of experimental data and 
the thermal production can be attributed to prompt and pre- 
equilibrium contributions. 



thickness function of the nucleus A; 

f +00 



dzpAix,y,z) , 



(3) 



where the nuclear density is given by a Woods-Saxon 
distribution, 



PAir) 



Po 



1 + exp [(r - R) /a] 



with the normalization. 



d?r pA{r) = A. 



(4) 



(5) 



A similar expression holds for the nucleus B. We shall 
take the nuclear radius R for the Pb nucleus to be 6.5 fm 
and the diffuseness a to be 0.54 fm. The nucleon-nucleon 
inelastic cross-section a is set to 32 mb relevant for the 
= 17.3 GeV. 

Now we assume that the effective number of protons 
from the projectile or the target which constitute the 
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FIG. 11: Spectra (upper pannel) and V2{pt) (lower panel) ' 
for primary pious from Pb + Pb collisions having 6 = 7 fm at 
SPS energy, for different initial times. 



participants at a given impact parameter b is given by; 

^^part(&) 



^Proj - ^Targ ^ ^ 



(6) 



with a similar expression for the effective number of neu- 
trons. These numbers then decide the shadowing func- 
tions Ra„ii{x, Q^) as well as the effective structure func- 
tions. 



fA^aix) ^ ^fpi.x) -f ^fnix) (7) 

^cff ^cff 



for the prompt photon calculations. We multiply the 
cross-sections for the production of photons from the 
"effective" nucleon-nucleon collisions using NLO pQCD, 
with the nuclear overlap function; 



Tab(6)- j dxdyTA{x+^-,y)TB{x-^-,v) (8) 



to get the yield of prompt photons. Finally the yield is 
averaged over the impact-parameter range covered by the 
centrality of the collision. 
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FIG. 12: Spectra (upper panel) and U2(pt) (lower panel) for 
primary rho-mesons from Ph + Pb collisions having & = 7 fm 
at SPS energy, for different initial times. 



FIG. 13: Spectra (upper panel) and V2{pt) (lower panel) for 
primary protons from Pb + Pb collisions having 6 = 7 fm at 
SPS energy, for different initial times. 



III. THERMAL PHOTONS 
A. Initial conditions 

We have already noted that the importance of ther- 
mal photons lies in their sensitivity to initial conditions. 
The simplest and most widely used initial conditions as- 
sume formation of a hot, thermalized, and chemically 
equilibrated quark gluon plasma at some initial time tq, 
beyond "which the system expands isentropically ignor- 
ing the viscosity effect. This makes the powerful meth- 
ods of hydrodynamics available to us. One may also use 
a parametrized fire-ball to describe the evolution of the 
system. 

For this study, we employ a boost invariant hydrody- 
namics dl] as our model of the underlying bulk evolution, 
especially for the purpose of obtaining the initial energy 
and temperature distributions. This model has been used 
extensively to explore and hadron pro duction and elliptic 
flow of hadrons as well as photons wi\ and dileptons z9| . 
For the SPS energies under consideration, the initial con- 
ditions are estimated by assuming [l^, [2^ [27|, [111 that the 
deposited energy in the transverse plane is proportional 
to the number of wounded-nucleons [29| (or participants); 

e{x,y,b,To) = Kv{x,y,b) (9) 



where is a constant and i^ix, y,b) is given in Eq. ^ . We 
further assume, as in Ref. [23] that the initial transverse 
density profile of net baryon number is proportional to 
the participant profile as well 

n{x,y,b,TQ) ^ Le{x,y,b,To) . (10) 

The authors of Ref. [1^ have shown that taking tq = 
0.8 fm/c along with K = 2.04 GeV/fm and L = 0.122 
GeV~^ provides a remarkably quantitative description of 
the particle spectra measured for the Pb+Pb collisions at 
V^NN ~ 17.3 GeV, using the equation of state Q, which 
provides that the thermally and chemically equilibrated 
QGP undergoes a first order phase transition to hadrons 
at Tc ~ 164 MeV. We have checked that these values 
give a quantitative description of the deposited transverse 
energy measured by the WA98 experiment [3^ for central 
collisions. 

In the present work we use these values for K, L, and 
Tq. We additionally explore the consequences of varying 
To such that the rapidity density of total entropy, dS/dy, 
and net baryons, dN^-g/dy, remains fixed (see Ref. |28l. 
|34| for a similar approach). This is attained by taking 
the entropy density oc e^^^ and then using sqTq and uoTq 
as constants, where Sq — s{x — 0,y — 0,b — 0) and 
uq = n{x = 0,y = 0,6 = 0). This corresponds to an 
isentropic expansion. Note that the shape of the these 
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distributions are taken as independent of tq. Thus we 
have used tq = 0.2, 0.4, 0.6, 0.8, and 1.0 fm/c. These 
then correspond to the peak temperature, To{x — 0,y — 
0, 5 = 0) of 420, 330, 284, 257, and 238 MeV respectively, 
while the average temperatures are 380, 295, 257, 233, 
and 215 MeV, respectively. 



B. The Flow Patterns 

As a first step we determine the time-evolution of the 
average energy density, average temperature, and the av- 
erage transverse velocity of the expanding system ob- 
tained from the hydrodynamic calculations for a central 
collision. We take the average by defining. 



(/)- 



/ dxdy f{x,y)e{x,y,T) 
J dxdye{x,y,T) 



(11) 



We note (see Fig. [3]) that in the over-lapping time-span, 
the variations of these quantities are quite similar, though 
an earlier start leads to a slightly larger build-up of the 
flow velocity and a faster cooling of the system. In the 
final stages the temperatures and the velocities do not 
differ beyond about 10% for different initial times, though 
the energy density varies by about 40% (as it varies as 

Photons are sensitive to the initial temperature. 
Therefore, an earlier initial time with higher initial tem- 
perature will lead to a considerably enhanced production 
of photons at higher transverse momenta. On the other 
hand, this should not affect the spectra of hadrons since 
they are emitted at much later freeze-out stage when the 
effect of having different initial times is mostly washed 
out. 

In Fig. [4] we have shown the time-evolution of the av- 
erage effective temperature (or the blue-shifted temper- 
ature) to see the combined effect of the cooling and ex- 
pansion (velocity). We define the effective temperature 
as, 



T 



1 -I- wr 

1 — Vt 



(12) 



We note that as in Fig. [31 the results differ only 
marginally beyond the time of about 1 fm/c, confirming 
once again our surmise that the difference in the produc- 
tion of thermal photons should mostly arise from contri- 
butions before this time. 



C. Thermal photons 

We calculate the production of thermal photons by 
folding the history of the evolution of the system with the 
rate for the production of photons from the quark matter 
and the hadronic matter. We use the complete leading- 
order results for the production of photons from the QGP 
from Arnold, Moore, and Yaffe i31il and the latest results 



for the radiation of photons from a hot hadronic gas ob- 
tained by Turbide, Rapp, and Gale |20|[. As mentioned 
earlier, the equation of state (EOS Q 2^) incorporat- 
ing a phase transition to quark gluon plasma at T « 164 
MeV, and resonance gas for the hadronic phase below the 
energy density of 0.45 GeY/im^^ is used to describe the 
evolution. The mixed phase is described using Maxwell's 
construction. The freeze-out is assumed to take place at 
e = 0.075 GeV/fm^. Final results are obtained by taking 
an average of the results over the range of impact pa- 
rameters b between and 4.6 fm corresponding to 0-10% 
most central collisions, considered by the WA98 experi- 
ment. 

We summarize our results for the case of tq =0.2, 0.4, 
0.6, 0.8, and 1.0 fm/c in Figs. l5l6l7l8l9l As expected from 
the discussion earlier, we find that the hadronic matter 
contribution to the single photons is only marginally al- 
tered as we increase the initial time or decrease the initial 
temperature. The quark matter contribution at large px 
however drops as the initial time is increased (the initial 
temperature increased) . 

We note that the prompt photon production is about 
17% of the total yield. Noting that these are NLO results, 
the lowest order prompt photon production is perhaps 
only of the order of 10% of the total single photon pro- 
duction measured in the experiment. We also note that 
the thermal production of photons is almost identical to 
the prompt photon production when tq — 0.4 fm/c. 

We have also shown the results for the " Thermal 
X Prompt" photon contribution, with k adjusted to re- 
produce the experimental results at pr = 2.55 GeV/c. 
It is good to see that the same normalization provides 
a good description to the entire pt range in every case 
(see middle panels). We find that scaling the prompt 
photon results by factors of 2.7, 4.9, 5.4, 5.7, and 5.9, re- 
spectively are necessary in order to provide a quantitative 
description of experimental results. 

We can perhaps argue that k accounts for the Cronin 
effect in case of nucleus-nucleus collisions as well as 
pre-equilibirum contributions which must surely be ac- 
counted for when tq is large. We do know that the pre- 
equilibrium electromagnetic radiations look thermal in 
nature [3^ . and we have noted that in the present case, 
the prompt and the thermal contributions have similar 
slopes, for large initial temperatures 

Even though a value of tq =0.2 fm/c may be consid- 
ered too small, let us not forget that in this notation the 
nuclei would start interacting at t = — i?/7 or at about 
T « —0.6 fm/c, and thus a hot and dense system can be 
considered to be formed soon after the complete overlap, 
which occurs at t = fm/c. 

Let us try to see if some additional experimental result 
could actually distinguish between the different values for 
To, and thus in a potentially interesting observation, we 
note (see lower panels of Figs. [51617181 [9)). that the elliptic 
fiow parameter V2 for the thermal photons [23 | is quite 
sensitive to the formation time tq [33| . We also note the 
peak at low pr in the V2{pt), first noted by authors of 
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Ref. (2J| and interpreted as a consequence of competi- 
tion of TTTT —^ p7 and np — > reactions. We note that 
as we decrease tq, the contribution of the quark matter 
increases. As this contribution arises from earher times, 
where the momentum anisotropics are smaher, decreas- 
ing To thus leads to an overaU reduction in V2 for ther- 
mal photons. Even though the azimuthal anisotropics 
for these fairly central collisions are small, they reveal an 
important sensitivity to the formation time (see Fig. [TOl 
Note also the inversion of order of the results for V2 with 
increasing tq with and without accounting for the non- 
thermal component). 

Let us pause here to consider a question which has 
troubled the analysis of single photons from the WA98 
experiment, from the very beginning. In the present 
work, we have started with the assumption of a forma- 
tion of QGP at time tq. However, several studies have 
also [13, [ss'l presented a reasonable description of the 
data by assuming only the formation of a hot hadronic 
gas in the collision without ever forming a QGP. Which is 
the right scenario, then? We note here that the photon V2 
provides a possible resolution to this question. If no QGP 
is formed, then the V2 for thermal photons will closely fol- 
low the V2{pt) for p mesons at larger px- Hence, it will 
be considerably lar ger than our prediction and also will 
rise monotonically [2J| as pr increases. This suggests 
that a measurement of the V2 of thermal photons along 
with their spectra could very firmly distinguish between 
the two scenarios. 

Coming back to our present discussion, we note that 
the results for V2 for direct photons will be modified from 
the values for the thermal photons due to the presence 
of prompt photons (see Fig. [TOl lower panel). However, 
the prompt photons as well as the pre- equilibrium pho- 
tons will not contribute to the azimuthal anisotropy of 
the photon distribution, as they are not subjected to any 
collectivity. We can safely neglect the small effect of az- 
imuthal dependence of jet-quenching which may affect 
the fragmentation photons (which is less than about 30% 
of the prompt contribution in the present case) or those of 
jet-induced photons [s^, which measure the anisotropy of 
the initial state [S^]. This is because the QGP, if formed 
at the SPS, is very short lived and not very hot, as in- 
dicated by a small jet-quenching (not exceeding about 
25-30%) for such collisions We show the results 

for final V2 for single photons for the case of the WA98 
experiment in Fig. 1101 

We finally recall a completely different calculation us- 
ing the parton cascade model [H, H^l for this case, where 
the scattering and radiating partons produced a not-so- 
dense partonic system, but it was enough to reproduce 
the single photon production seen by the WA98 experi- 
ment beyond about 3 GeV/c, if the partonic distributions 
are given an intrinsic < A:t > of about 0.44 GeV/c. In 
absence of the intrinsic < kx > the production is smaller 
by a factor of about 2. These photons can be considered 
as due to prompt and the pre-equilibrium contributions. 

We have so far assumed that the prompt and pre- 



equilibrium photons and their enhanced production due 
to intrinsic kx can be estimated by using a multiplica- 
tive factor K to the NLO pQCD results. Within this ap- 
proach, we have ascertained that the photon observables, 
especially the V2 , are sensitive measures of the initial con- 
dition. Admittedly, there are some uncertainties in our 
approach such as the NLO contribution to the thermal 
photon production and the effect of the viscous hydrody- 
namic evolution. At present, their effects are unknown 
although the effect of the finite viscosity on photons may 
soon be calculated [4l| . 



D. Particle spectra 

How will the reported good description of particle spec- 
tra [11] obtained using tq = 0.8 fm/c, be affected, if 
a different value is used for tq? Instead of discussing 
a complete calculation (with resonance decay accounted 
for), we show the primary spectra of pions, rho mesons, 
and protons for a typical impact parameter b =7 fm, for 
different values of tq, but keeping the entropy fixed as in 
the calculations discussed above (see Figs. lllll2ll3p . We 
note that as the inverse slope for all the cases rises with 
decrease in tq as the radial fiow sets in earlier. We have 
checked that the increase in the inverse slope for pions 
is about 11%, about 15% rho mesons and protons as the 
initial time is decreased from 1 fm/c to 0.2 fm/c. We 
also note that the change in the spectra for the primary 
particles is quite marginal for pT below 1.5 GeV/c, even 
though it varies by a factor of about 2> aX px — ?> GeV/c. 
What is most interesting is that the differential elliptic 
ffow parameter V2 for hadrons is almost independent of 
the initial time. 

We conclude then, that the good description of 
hadronic spectra for low transverse momenta at SPS en- 
ergies will remain unaffected by the reduction of initial 
time from 1 fm/c to 0.2 fm/c. This is in contrast to what 
we saw earlier for thermal photons. 



IV. SUMMARY AND CONCLUSIONS 

We have re-analyzed the single photon production in 
Pb + Pb cofiisions at the CERN SPS energies for 10% 
most central collisions. Several improvements have been 
incorporated. The iso-spin, shadowing, and impact pa- 
rameter dependence of the prompt photon production 
are explicitly included. NLO pQCD calculations are per- 
formed with the factorization, fragmentation, and renor- 
malization scales fixed at pt/2 based on a global descrip- 
tion of the available data for pp collisions. For the ther- 
mal photons calculations the initial conditions are taken 
as those which provided a good description to hadronic 
spectra, with a tq = 0.8 fm/c. We explored the conse- 
quences of using smaller initial times, keeping the entropy 
and the net-baryon number fixed. 
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We find that the data can be explained using a small 
formation tq of the order of 0.2 fm/c when supplemented 
with prompt photons evaluated at NLO pQCD with a 
K factor « 2.7 to account for the Cronin effect. Larger 
initial times require much larger values for k, which may 
be mimicking the pre-equilibrium contribution. A unique 
sensitivity to the formation time is seen in the photon el- 
liptic flow, which could be useful in ascertaining whether 
a QGP was formed at the SPS energy. 
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